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Abstract 

Motivated by the Ginzburg-Landau theory of superconductivity, we estimate in the semi-classical 
limit the ground state energy of a magnetic Schrodinger operator with De Gennes boundary con- 
dition and we study the localization of the ground states. We exhibit cases when the De Gennes 
boundary condition has strong effects on this localization. 
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I. INTRODUCTION 



Let J] C I 2 be an open bounded domain with regular boundary. Let us consider a 
cylindrical s_ d uc t i„ g sample of cross section «. The S „perconduc t i„ g proxies are 
described by the minimizers ("0)^4) of the Ginzburg-Landau functional (cf. Refs. 



)OC 

B 

J ||(V - lanA)^ 2 + crV|curlA - 1| 2 + y(|V| 2 - 1) 2 | dx 

+ [ 7|V(x)| 2 d/i| af2 (x), (LI) 
Jen 



which is defined for pairs (ip, A) G C) x M 2 ). The parameter k is a characteristic 

of the material. A material is said to be of type I if k is sufficiently small and it is said 
to be of type II when k is large. The parameter a is the intensity of the applied magnetic 
field which is supposed to be constant and perpendicular to Q. For a minimizer (if), A) of 
the energy Q, the function i/j is called the order parameter and \ip\ 2 measures the density of 
superconducting Cooper electron pairs; the vector field A is called the magnetic potential 
and cmLA is the induced magnetic field. Note that the order parameter ip satisfies the 
following boundary condition proposed by De Gennes^ : 

v ■ (V - ianA)^ + 77/; = 0, (1.2) 

where v is the unit outward normal of dQ and 7 6 R is called in the physical literature the 
DeGennes parameter. Note that the boundary condition (jl.2|) was initially introduced in 
the theory of PDE by Robin. 

The physicist De Gennes^ introduced the parameter 7 in order to model interfaces between 
superconductors and normal materials. In that context, 7 is taken to be a non-zero positive 
constant and - (called the extrapolation length) usually measures the penetration of the 
superconducting Cooper electron pairs in the normal material. The size of 7 depends on 
the nature of the material adjacent to the superconductor and it ranges from 7 = (inter- 
faces with insulators) to 7 = +00 (interfaces with magnetic and ferromagnetic materials). 
Experiments show that for superconductors adjacent to ferromagnetic materials, the order 
parameter if) vanishes at the boundary^ and the boundary condition (jl.2|) is changed to the 



Dirichlet boundary condition. Negative values of 7 were also considered in the physical 
literature^. It is suggested that negative values of 7 would be useful for modeling the sit- 
uation when a superconductor is adjacent to another superconductor of higher transition 
temperature. 

Suppose that we have a type II superconductor (i.e. k is large). The functional Q has a 
critical point of the type (0, A). Such a critical point is called a normal state. It is then 
natural to study whether a normal state is a local minimum of Q in the presence of a strong 
applied magnetic field. The Hessian of Q near a normal state is given by : 



V - ianA)(f)\ 2 - K 2 \(f)\ 2 ) dx + / 7 \(j)\ 2 d^\ dn {x) + (an) 2 / |curl.B| 2 cfx 

Jan Jn 



By defining the change of parameter h = — , we have then to study as h — > the positivity 
of the quadratic form : 



H\Q) 3u^ II (W - iA)u\\ 2 L2{n) + h 2 7 \u\ 2 d^\ an {x) - {k}i) 2 \\u\ 

Jan 
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The semi-classical limit h —>■ is now equivalent to a large field limit a —>■ +00. In order to 
study the influence of the size of 7, it seems reasonable to suppose that 7 is depending on h. 
Also, due to the possibility of having different materials exterior to Q together with possible 
lack of symmetry in the geometry of Q, it seems also convenient to take 7 as a function 
of the boundary. Thus, given a vector field A G C°°(fi;IR 2 ), a regular real valued function 
7 G C°°(dfl] R) and a number a > 0, let us define the quadratic form : 

B u _> q ^ {u) = \\ {hV _ iA)u\\ 2 L , {n) + h 1+a f 7(*)K*)| 2 % 5 >)- (1-3) 

Jan 

Observing that n is semi-bounded, we consider the self-adjoint operator associated to 
Ih A n by Friedrich's theorem. This is the magnetic Schrodinger operator P^'J n with domain 
D(P$ fl ) defined by : 

KXn = -{hV-iA) 2 , 

D(P^ln) = {u^H 2 (n); v ■ (W - iA)u lan + h^u ]gn = 0}. 
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We denote by /^(a, 7, /i) the ground state energy of P^'Jq which is defined using the 
min-max principle by : 

a, 7 / \ 

(1)/ n • r 1h,A,n\ U ) /j r\ 

/i l j (a,7, h) := mf —-^ . (1.5) 

ueffi(n),«^o ||«||2a( n j 

Let us recall also that this eigenvalue problem is gauge invariant. 

In the case when 7 = (which corresponds to a superconductor surrounded by the 

vacuum) , a lot of papers are devoted to the estimate in a semiclassical regime of the ground 

state energy of P^aq- We would like here to mention the works of Baumann-Phillips-Tangfi, 

Bernoff- Sternberg^, del Pino- Felmer-Sternberg 8 , Helffer-Mohamed^, Helffer- Moraine^ and 

the recent work of Fournais-Helffer 11 . The special case when a = 1 and 7 is a positive 

constant was considered by Lu-Pa n 12 i 13 . It was shown that in this case the effect of the 

u (1) (l,7, h) 

De Gennes parameter 7 is weak in the sense that the limit lim ^ — - — is the same as in 

h^O h 

the case 7 = 0. This regime is therefore not sufficient to recover all the physically interesting 
cases considered in Refs. 3 U It 

in Ref. Q and to analyze (for all values of a) the influence of the boundary term in (|I.3J) on 
the localization of the ground state energy of the operator P^'Jq- 

Following the technique of Helffer-Morame^, we have to understand the model case of 
the half-plane when the magnetic field and the function 7 are both constant. Consider the 
magnetic potential : 

A (x 1 ,x 2 ) = « ( _a; 2, xi), V(ari, x 2 ) G R x R + . (1.6) 
Notice that curlA = 1. Let us define the function 

R 3 7 ^ 6(7), 

where 



is the object of this paper to establish the results announced 



v - ^oHI! 2(KxR+) + 7 Jr K27, o)| 2 ob 



and 



6(7) := inf , ,,7 J (1.7) 

U effi o (RxR + WO \\ U \\h(RxU + ) 

H\ o (R x R + ) = {u G L 2 (R x R + ); (V — iA )u G L 2 (R x R + )}. (1.8) 



Note that 0(7) is the bottom of the spectrum of the operator P^'J q with h = 1 and 
= Rx R + . We shall see that ©(7) < 1 (cf. Theorem ITL2J) . If 7 = 0, we write : 

6 := 6(0). (1.9) 

It is Go which appears in the analysis for the Neumann proble m 6 ' 7 ! 8 ! 9 ! 10 ! 11 ! 12 ' 13 . Actually, 
we are interested in the bottom of the spectrum of the operator -P/^a ,rxir + but a scaling 
argument gives us : 

Vh G R + , Vq, 7 G R, inf Sp (P^ ,RxR + ) = ^6(^ 1/2 7). (LIO) 

The semiclassical analysis of the half-plane model depends then on the sign of both a — ~ 
and 7. We have then to investigate the asymptotic behavior of 0(7) when 7 — > and when 
7 — > ±00. This will be the object of study in Section ITT1 

Now we state our main results. 

Theorem 1.1 Suppose that Q C M 2 is open, bounded, connected and having a smooth bound- 
ary. Suppose moreover that the magnetic field is constant curl A = 1. Then, for a > and 
7 G C°°(<9fi;R), the ground state energy of the operator P£'J n satisfies : 

/x«(a, 7, /*) = he(h Q ^ 2 l0 ) (1 + o(l)) , (h - 0), (1.11) 

w/iere 70 := min xe9n 7(2;). 

Theorem 11.11 gives a first term approximation of fj,^(a,j,h). The asymptotics (jl.ll|) is 
valid without the need to any non- degeneracy hypothesis on the set of minima of 7, and 
holds for the function 7 being constant as well. Let us remark that the asymptotics (jllljl 
depends strongly on a. In particular, when a = |, we get : 

I™, r = © 7o < 1, 

and if 70 = or if a > |, then (cf. Proposition III.5|) : 

7, /i) 

hm = B < 1. 

h-»o ft, 
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When a < |, it is the sign of 70 that affects the asymptotics. Actually, if 70 < we have 
(cf. Proposition 111.8)1 . 

/i (1) (a,7,/i) 2 
hm — = -7 , 

h^o h 2a 



and if 70 > 0, we have (cf. ()II.46)) ) 



/x (1) (a,7,/i) = 



lim 

h->0 a 

which is the same behavior as that for the Dirichlet problem^. This 



ast regime (0 < a < \ 



and 70 > 0) is in accordance with the physical observations in Ref. ]4 . 



In the next theorem, we give a two-term asymptotics of fi^(a, 7, /i) when a g]^, 1[. 

Theorem 1.2 Suppose in addition to the hypotheses of Theorem \I. 1\ that 5 < a < 1 and 
t/iat t/ie function 7 zs non- constant. Then we have the following asymptotic expansion as h 
tends to : 

f]W(a,i,h) = hQ + 6M 37o /i a+1 / 2 + C(/i inf ( 3 / 2 ' 2a )), (1.12) 
where M 3 a strictly positive universal constant. 

The constant M3 satisfies O'(0) = 6M3 and it will be defined precisely in Section |Hj 



see however fill. 25)) and pi. 27)) . Comparing with the result obtained in Ref . [l0|, the second 
term in the two-term asymptotics of ^'(a, 7, h) when 7 = is of order /i 3//2 , whereas it 
is of order h a+1 l 2 in the regime considered in Theorem 11.21 Let us mention also that in 



Ref. 



11 



the authors obtain (when 7 = 0) a complete asymptotic expansion under a generic 
hypothesis on the scalar curvature of dVt. It seems that a complete asymptotic expansion 
could be obtained in the regime of Theorem 11.21 but under the following generic hypothesis 
over 7 : 

- 7 has a finite number of minima; 

- all the minima of 7 are non-degenerate. 

We leave this point hoping to analyze it in a future work. 
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Next we turn to the question of the localization of the ground states. Let be a 

ground state of the operator P^ n . We say that u a ,y,h is exponentially localized as h tends 
to near a closed set B in Q if there exists f3 > 0, and for each neighborhood V of £>, there 
exist positive constants h , S and C such that : 



\U a ,-y,h\\L 2 (n\V) 



< 



Cexp (^~^pj ||«a, 7 ,fclU 2 (fi), V/ie]0,/io]. (1-13) 



In the next theorem we describe some effect of 7 on the localization of the ground states 
of the operator P^'Jq- 

Theorem 1.3 Under the hypotheses of Theorem \I.l[ if j < or | < a < 1, a ground state 
of the operator P^J n is exponentially localized as h tends to near the boundary points 
where 7 is minimum. 

More precisely, \1.13\) is satisfied with [3=1 — a if jo < 0, j3 = if | < a < 1, and 
[3 = 1/2 otherwise. 

In the special case a = 1, the scalar curvature k t and the function 7 affects the asymptotic 
expansion of the ground state energy to the same order. 

Theorem 1.4 Suppose in addition to the hypotheses of Theorem \l.l\ that a = 1. Then we 
have the following asymptotic expansion as h tends to : 

/i«(a, 7 , h) = hQ - 2M 3 ( Kr - 3 7 ) m ax/> 3/2 + 0(/* 13/8 ), (1.14) 

and a ground state u an ,h of the operator P h '^Q is localized near the boundary points where 

the function k t — 37 is maximal. 

More precisely, U.l^l is satisfied with [3 = 1/4. 

If 7 is constant, the remainder in (jI.14j) is better and of order 0(h 5 ^ 3 ). When 7 = 
we recover in the above theorem the result of Helffer-Morame^. Let us mention that the 
expansion (II14J) is announced by Pan 25 in the particular case when 7 is a positive constant. 



As in Ref. 



11 



we believe that an asymptotic expansion with higher terms could be obtained 
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under a generic hypothesis on the function k v — 37. 



In the next theorem, we study the case when the function 7 is constant and we find that 
only the scalar curvature plays a role. 



Theorem 1.5 Suppose in addition to the hypotheses of Theorem that the function 7 
is constant and that a > |. There exists a constant M 3 (a;, 7) > such that we have the 
following asymptotic expansion as h tends to : 

M (1) («,7, h) = /i9(/i a - 1 /2 7 ) - 2M 3 (a, 7 )( S .) max /i 3/2 + o(h^ 2 ). (1.15) 

Moreover, a ground state of the operator P^'J n is localized as h tends to near the boundary 
points where the scalar curvature is maximal, and HI. 1 3\) is satisfied with (3 = 1/4. 



a = i 



When a > |, the constant M 3 (a,7) is equal to the universal constant M 3 . When 
we have M 3 (~,7) = M 3 (7), where the constant M 3 (7) will be defined in Section II (cf. 

(HE!). 

This paper is organized in the following way. In Section |Hj we link the analysis of the 
half-plane model operator to that of a one dimensional operator. We get in particular 
the existence of a number £(7) > such that 0( 7 ) is the lowest eigenvalue of the operator 
— <9 2 + (t — £( 7 )) 2 . Let <p 7 be an eigenfunction associated to 6(7). We establish the regularity 
of 0(7) and ip 1 as functions of 7, the asymptotic behavior of 6(7) as 7 — > ±00, and uniform 
estimates with respect to 7 describing the exponential decay of tp 1 at infinity. 



In Section ITTT1 we use the eigenfunction </? 7 to construct a test function inspired by Refs. 



10 



and we obtain an upper bound for ^(a, 7, h). We then carry out a similar analysis to that 



in Ref . [10j and we use the results of Section UH to prove Theorem 11.11 
In Section IIV1 we show how to get the localization of the ground states using Agmon's 
technique^. Finally, in Section the analysis of a one-dimensional family of operators on 
a weighted L 2 -space appears (cf. (|V.21J) ). It is the same family of operators appearing in 



Ref. 



lOj (Section 11) but with a different boundary condition this time. This analysis permits 



S 



us to derive two-term asymptotics of the ground state energy showing the influence of the 
scalar curvature. We finish then the proofs of Theorems 11.21 1131 11.41 and 11.51 

II. THE MODEL OPERATOR 

Given 7 e R, let us consider the quadratic form : 

H 1 Ao {RxR + )3u^q[ 1 ]{u) = \\{V-iA )uf LHM ) + 1 [ Ka^O)! 2 ^!. (II. 1) 

Jr 

The magnetic potential A and the form domain H\ (M. x R + ) are defined respectively 
in (jl.6|) and ()L8[) . Observing that the quadratic form q[y] is bounded from below, we can 
associate to q[j], by taking the Friedrichs extension, a unique self-adjoint operator P[j] on 
L 2 (R x R+). The min-max principle gives that the bottom of the spectrum of P[j] is equal 

to 0(7) (cf. jTTD). 



A. Link with a one dimensional operator 

By a change of gauge and a partial Fourier transformation with respect to the first 
variable, we obtain that the spectral analysis of the operator P[j] will be deduced from that 
of the ^-family of one dimensional operators : 

^[7,e] = -^ + (^-0 2 , (H.2) 

with domain 

D(H[y,£]) = {u G £? 2 (R+); u'(0) = 7 «(0)}, (II.3) 
where, for a given integer fc, the space B k (R + ) is defined by : 

B k (R+) = {ue H k (R + ); t k u G L 2 (M + )}. (II.4) 

Note that the operator H [7, £] has compact resolvent and hence the spectrum is discrete. 
We denote by ^^(7, C) the first eigenvalue of ifpy, £]. The min-max principle gives : 

ll m M= taf 



where g[ 7 ,£] is the quadratic form associated to H[j, £] : 



£](«)= / {\u'(t)\ 2 + \(t-Ou(t)\ 2 )dt + 1 \u(0)\ 2 . 
A spectral analysis using the separation of variables (cf. Ref. \lw gives us 

e( 7 ) = mf/i«( 7 ,e). 

In the following lemma, we collect some useful estimates of //W(7, £). 



Lemma II. 1 Given e g]0, 1[, we /iave, 

^ (1) (7,0>(i-^ (1) (o,0 

t 

where 7 _ = max(— 7 , 0). 
Moreover, given 7 G IR ; we /iai>e : 

lim /i (1) ( 7 ,0 = +oo, lim /i (1) ( 7 ,0 = 1- 
4— >— 00 £^+00 

Proof. Using the density of C7(l~) in H 1 (R+), we get for any u G H 1 (R + ) 

POO 

\u(0)\ 2 = -2 u(t)u\t)dt. 



By the Cauchy-Schwarz inequality, we get for any a > : 

|u(0)| 2 < a\\u\\ 2 L2{R+) + -|K||i 2(R+) . 
Taking a = - (with 7 < 0), we get : 

? [ 7 ,£](«)>(l-e)g[0,£](« 

The min-max principle now gives (|II.7j) . 



(7-; 



Notice thai (ICS) is valid for - (Ref.Q). So the limit as £ — > — 00 m i|ii.g|) is now a 
consequence of the estimate (|II.7|) . For the reader's convenience, let us give for non-zero 7 
a proof for the limit as £ — > +00 in (|II.8[) . Let us denote by /i D (£) the first eigenvalue of 
the Dirichlet realization of the harmonic oscillator — d 2 + (t — £) 2 on M + . We have by the 
min-max principle : 



(II.5) 



(II.6) 



(II.7) 



(II.8) 



(II.9) 



(11.10) 



^ (1) (0,O +7l^, f (0)| 2 < ^ {1) (7,0 < /AO, 



(11.11) 
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where cp y ^ is the L 2 -normalized eigenfunction associated to fi^ (7 , £) . Let us notice also 
tha1^i& lim n D (£) = 1. So, if we know that lim |y? 7l? (0)| 2 = 0, then (jTTT"T|l is sufficient 



to deduce the limit as £ — > +00 in ()II.8jl . Thus, it remains for us to prove the following 
claim : 

Given e GlO, 1[ and 7 G M, there exists a constant C > such that, 

(11.12) 

V£e[C,+oo[, |^, € (0)| 2 < Ce"^/ 2 . 

Let us mention that the decay in (j!I.12|) is not optima l 17 ' 18 . We prove (jII.12|) using Agmon 
type estimated. Let $ be a regular function with compact support. An integration by 
parts gives the following identity : 

g[ 7) £] (eV?) = ^(7,0 ||eVdl^( R+ ) + ll $ ' e Vdl^ +) ■ ( IL13 ) 

Using the estimate (jll.lOj) (with e = 1/2) together with the fact that /i D (£) is bounded for 
£ G R+, we can rewrite (|II.13j) in the form : 

\ (l ( e *^)'|| L2(R+) + IK* " - ^ ll e *^||l 2(R+) + ||*'e*^^||| 2(M+) , 

(11.14) 

for some constant C > 0. We choose now $ as : 



$(t) 



^(i^l! i/ < t < 1, 
if t>l. 



Under this choice of $, we can get a sufficiently large constant C > such that, for £ G 
[C, +00 [, we can rewrite ()H.14|) in the form : 

ll e *^7,dli/i(n) - c - 

Using the Sobolev imbedding if 1 (M + ) ■— > L°°(M + ), this last estimate is sufficient to deduce 

(irmi) . □ 

Following the analysis of Dauge-Helffei^, we have now the following result. 
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Theorem II. 2 For each 7 6 I, 6(7) < 1 and the function I9(h /iW(7,£) attains its 
minimum at a unique positive point £ (7) that satisfies : 



£(7) 2 = e( 7 ) + 7 2 



n 



(11.15) 



Proof. Let us notice that by Kato's theory (Ref. 120), the maps 



£^«( 7 ,£), («^^ 2 W 

are analytic. Here we recall that tp^^ is the unique strictly positive and L 2 -normalized 
eigenfunction associated to ^'("f, £). Let us consider r > 0. Note that, 

/i« (7, £ + r)<p l4+T (t + r) = If [7, £] (<f l4+T (t + r)) , W G R + . 

Taking the scalar product with ip y ^ and then integrating by parts, we get : 



(11.16) 



( yU «( 7 ,£ + r)-^ 1 )( 7 ,0) / ^ +T {t + r)^{t)dt 
= ^, £+t (t)^(0) - 7^ 7i5+T (r)^ 7i5 (0). 



Recall that we have the boundary conditions : 



^,£+r(0) = 7<^7,C+t(0), ¥\e(0) = 7^(°)- 



Then we can rewrite f|11.16|) as : 

/i (1) (7,e + r)-/iW(7,0 



(p^ +T (t + r)ip^(t) dt 



^(0). 



By taking the limit as r — > 0, we get : 



^ (1) (7,O = « f (0)-7< e (0))^(0). 



Finally, we make the substitutions : 



< c (0) = (e 2 -A* (1) (7, 0)^(0). ^, e (0) =7^(0), 
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and we get the following formula, 

^«( 7 ,e) = (e 2 -^ (1) (7,e)-7 2 )i^(o)i 2 , (im 



called usually the F-formula (cf. Refs. Il9lj2ll ). Using (jll.7|) and (jll.8|) . we get : 



%/ (1) ( 7 ,O|^o<0, (7, >0, 

for a sufficiently large r] > 0. This gives the existence of a positive critical point of ^'('j, £). 
Let us notice now that for any critical point £ c of /x^(7, £), we have : 

df//«( 7 ,£)|^ c = 2£ c |y^(0)| 2 . 

This shows that any negative critical point is a global maximum and any positive critical 
point is a global minimum of /x^(7, £)• Coming back to (jll.8|) . lim //^(7 , = +00, 
and thus there does not exist any negative critical points. Therefore, the minimum of 
£ 1— > /i^(7, is attained at a unique point £(7) > and the function £ 1— > /i < - 1 ^(7,£) is 
strictly increasing on [£(7),+oo[. This proves in particular (recalling (jll.6Jl ) : 

0(7) = ^ (1) (7,e(7))<l- 

□ 

In the sequel, we denote by y9 7 the unique strictly positive and L 2 -normalized eigenfunc- 
tion associated to the eigenvalue 6(7), and by H[y] the operator if [7,^(7)] : 

¥? 7 = ¥^ (7) , H[ 1 ] = H[ 1 ,ai)}- (H.18) 
In the next lemma, we collect various useful relations satisfied by the eigenfunction ip T 



These relations are similar to those given in Appendix A of Ref. 
Lemma II. 3 For each 76M, the following relations hold : 
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(t-£(7))M*)| 2 ^ = 0, (11.19) 
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(* - m?\^(t)\ 2 dt = - ||^ 7 (0)| 2 , (11.20) 



(* - m?WM 2 dt = \[1- 2(7£(7)) 2 ] l^ 7 (0)| 2 . (IL21) 

Proof. We follow the calculations done in Bernoff- Sternberg 7 . Let us consider the differen- 
tial operator : 

L = -d 2 t +{t-i( 1 )f-Q( 1 ). 
Note that for any polynomial p, we have the following identity : 

L(2 m - pV 7 ) = (P (3) - 4 [(t - £( 7 ) - 8(7)] V - 4(f - £( 7 ))p) <p T (11.22) 
Let v = 2p<£> 7 — p'ip-y. Integrating by parts we obtain : 

P+OO 

/ ^(t)(Lv)(t) dt = (v'(0) - 7^(0) V 7 (0). (11.23) 
Jo 

Taking p — 1, we get : 

/•+oo 

-4 / (f - e( 7 ))|^(t)| 2 rft = 2 (£( 7 ) 2 - 7 2 - 6( 7 )) K(0)| 2 . 



Recalling ()II.15|) . the above formula proves (|II.19|) . 

We prove (|II.20|) by taking p = (t — £(7)). To prove (jll.21j) . we take p= (t- £(7)) 2 . Note 
that we have in this case : 

t/(0)-7^(0) = 2(2( 7 £(7)) 2 -l) <p 7 (0). 

We get now from (ITT221 and (ITL231 : 

p+00 

-12 / (t - ai)?\^(t)\ 2 dt = 2 (2( 7 £( 7 )) 2 - 1) K(0)| 2 . 
Jo 

This proves (|iL2T|l . □ 

For 7 G R, let us define the parameter : 

M 3 (7) = ^(l + (7£(7)) 2 )l^(0)| 2 , (IL24) 
14 



and when 7 = 0, we write M 3 := M 3 (0). Note that ()11.21jl gives : 

M 3 = [ (t-£ ) 3 \<p (t)\ 2 dt, (11.25) 



where £0 := £(0). The constant M3 is the universal constant appearing in Theorems IL2 
II. 4[ and the parameter M 3 (7) appears as Ms(~,j) in Theorem 11.51 



B. Regularity 

We discuss now the reg u,a rity of the function. 7 - 9( 7 ) e > and 7 ~ , 7 6 It 
seems for us that Kato's theory (cf. Ref. 20) do not apply in this context at least for the 
reason that we do not know a priori whether the expression of the operator 

depends analytically on 7. Inspired by Bonnaillie 22 , we use a modification of Grushin's 
method 23 and we get the following proposition. 

Proposition II. 4 The functions R 3 7 h-> 0(7) G E and R 3 7 h-> <^ 7 G L 2 (R + ) are C 00 . 
Moreover, the function R 3 7 1— > <y9 7 G L°°(]ri + ) zs locally Lipschitz. 

The specific difficulty in proving Proposition 111.41 comes from the fact that both the 
expression and the domain of the operator H[y] depend on 7. To work with an operator 
with a fixed domain, we consider a cut-off x that is equal to 1 on [0, 1] and we apply 
the invertible transformation cp 1— » <p = e~ 7 * x ^V that transforms the boundary condition 
<p'{0) = 7^(0) to the usual Neumann boundary condition <p'(0) = and leaves the spectrum 
invariant (cf. Proof of Proposition 111.7)) . 

In the next proposition we determine 0'(7). This is a first step in the proof of Proposi- 
tion HH 
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Proposition II. 5 The function 7 1— > 6(7) is of class C 1 and satisfies : 



0'(7) = I^(O)| 2 . 



(11.26) 



In particular, we have : 



6'(0) = 6M 3 . 



(11.27) 



Remark II. 6 Using Formula ill. 15)) we get a/so £/ja£ the function 7 1— > £(7) zs 0/ c/ass C 1 . 

Proof of Proposition 111. 51 Let r be a real number. We shall define the following trial 
function : 



By standard Fredholm theory, the operator (#[7] — @(7)) _1 is defined on the orthogonal 
space of y? 7 and has values in D(H[j]). Hence, the function u\ is well defined, thanks to 
()11.19|) . and the function u satisfies the boundary condition u'(0) = (7 + r)u(O). When r is 
sufficiently small, it is a result of the exponential decay of y? 7 at +00 (cf. Propositions 111.91 
and 111.10)) and standard elliptic estimates that u 6 5 2 (IR + ). Therefore, u G D(H[~f + r]), 
and we have : 



u = e r *(<£> 7 + rui), 



where 



wi = - e(i)r l {^ 7 (o)|V 7 + V 7 + 2 (fa + r ) - *(-r))(* - fa))v,} • 



#[ 7 + r ]( u ) = e Tt {-d 2 t +{t- fa + t)) 2 - 2r<9 t - r 2 ) fa + ru x ). 



(11.28) 



Using the decomposition : 



[7 + r] = if [7] - 2(£( 7 + r) - £( 7 ))(f - fa)) + (fa + r) - fa)) 2 , 



we can rewrite (|11.28|) as : 



(#[ 7 + r ]_e (7 )_l^(o)|V) (u) 

= rV< (0( 7 ) + |^ 7 (0)| 2 ) Ml + (£(7 + r) - fa)) 2 u. 



(11.29) 
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We make the following claim : 

V 7 GM, 3C>0, We [-1,1], |e(7 + r)-e(7)l<C|r|. (11.30) 

Therefore, thanks to (|II.29jl and (jll.30|) . there exist constants C, r > such that, for all 
r E [—t , r ], we have : 

|| (#[ 7 + T ] _ 9 ( 7 ) _ |^ 7 (0)|V) u|| L2(H+) < CV 2 ||u|| L2(R+) . 

We get now by the spectral theorem the existence of an eigenvalue 0(7 + r) of the operator 
H[y + t] that satisfies the following estimate : 

|e(7 + r)-0( 7 )-|^ 7 (O)| 2 r| < CV 2 , Vre[-r ,r ]. (11.31) 

We make now another claim : 

V 7 GM, 3d>0, Vre[-l,l], |/i (2) (7 + r,e(7 + r))-/i (2) (7,^(7))| < C\t\, (11.32) 

where for (77, £) 6lxR, /i( 2 )(r/, £) denotes the second eigenvalue of the operator H[r},£\. 
Under the above claim, the estimate (|II.31j) gives : 

6(7 + r) = 6(7 + r), Vre[-r ,r ]. 

Consequently, we get that 6(7) is differentiable and satisfies formula (|II.26|) . We make now 
a final claim : 

The function 7 1— > |<£> 7 (0)| 2 is locally Lipschitz. (11.33) 
To achieve the proof of the theorem, we only need to prove (jll.30|) . (jll.32j) and (jll.33|) . 

Proof of num . 

As we have the formula (|H.15|) . it is sufficient to prove : 

V7GI, 3C>0, Vre[-l,l], |6(7 + r) -6( 7 )| < C\r\. (11.34) 

The min-max principle gives : 

^(7,0 + r|^ 7+Ti€ (0)| 2 < /x (1) (7 + r,0 < ^(7,0 + rK e (0)| 2 , V£ E R. (11.35) 
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Thus, given an eigenfunction ip of #[7, £], we need to estimate |</?(0)| 2 . Let u G D(H[j]). 
Using (jll.9|) we get : 

|u(0)| 2 < 2||n|| i 2 (R+) ||M / || i2(R+) . (11.36) 
We use now fllTTOl) (with e = 1/2) to obtain : 

\W'\\h(m + ) < 2g[7,e](«) + (7-) 2 |klli 2 ( K+ ). (H.37) 
Combining (|II.36|) and ()II.37|) . we get after an integration by parts, 

K0)| 2 < 2 ||#[7,£HL 2(K+) ||«|Ua (K+) + (l-) 2 \H\l HR+) , Vu G D(H[ 7 ]). (11.38) 
Let M := sup £(7 + r). Let us show that M < +00. Actually, the min-max principle 

re[-l,l] 

gives : 

^ {1) (7-l,0</^ {1) (7 + r,0<^ (1) (7 + l,0, Vre[-l,l], V(Gl. 
Recalling (jll.6|) . we obtain 6(7 — 1) < sup 6(7 + r) < 6(7 + 1). Formula (|II.15|) now 

re[-l,l] 

gives M < +00. Therefore, (jTT38j) gives : 

K, 5 (0)| 2 <C |^ 7+r , c (0)| 2 <C, Vre[-l,l], V£g[-M,M], 
for some constant C > 0. Consequently ()IL35|) yields the estimate : 

- CV_ < /i«( 7 + r,0 < /i (1) ( 7 ,0 + CV+, V£ G [— M, M). 
Minimizing with respect to £, we get (HHSII), thanks to Theorem HO 

Proof of num . 

Let u G S 1 (M + ). We shall compare q[j + r, £(7 + r)](tt) and q[j, £(7)](w)- In fact, we have : 

g[7 + r,£(7 + r)]M = g[7,£(7)]0) - 2(£( T + r) - £(7)) / (t - Z(i))\u(t)\ 2 dt 

Jo 

r+00 

+ (£(7 + r)-£( 7 )) 2 / \u(t)\ 2 dt + r\u(0)\ 2 , 
Jo 

where, combining ()II.36|) and (jll.37|) . 

K0)| 2 <2g[7,£]( M ) + (7-) 2 || M ||i 2(R+) . 
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hoo 

< / \( t -ai)Ht)\ 2 dt+ I \ u (t)\ 2 dt. 

'o Jo 



The Cauchy-Schwarz inequality gives : 

-+oo 

(t - m)Ht)\ 2 dt 

Using (jll.30|) . we get a constant C > such that, for all r G [—1, 1], we have : 

(l-CT_)g[ 7 ,e(7)]H-^_||«|| 2 i2(K+) < ff[ 7 ,£(7 + T)](u) 

< (1 + CY + )g[ 7 , £( 7 )](it) + Cr + \\u\\ 2 L2{R+) . 

The min-max principle proves now the claim. 

Proof of num . 

Let u = (p*y(t) — e~ Tt (p 7+T (t). It is sufficient to prove that : 

\u(0)\ 2 < C\t\, \/re[-r ,r }, (11.39) 

for constants C,r > 0. Using (|II.38|) . we have to estimate |M|l 2 (r + ) and ||-H"[7> £] u IU 2 (ir + )- 
Let / = (H[i\ - 6(7)) it. Then : 

/ = (e( 7 )-e( 7 + r))u + «;, 

where 

w(t) = e- Tt {-2rd t + 2(^(7 + r) - £( 7 ))(f - £( 7 )) - (£( 7 + r) - £( 7 )) 2 + r 2 ) <p 7 (t). 
Therefore, thanks to (jll.30|) and (jll.34|) . we have : 

\lHr + )<C\t\, y T e[-r ,r }. (11.40) 



Noticing that, after an integration by parts, (/, ^P 7 )l 2 (r+) — 0, we write, 

u = (H[ 1 ]-e( 1 )y 1 f. 

It is a standard result that the operator norm of (H[y] — 0( 7 )) _1 is bounded on the 
orthogonal space of </? 7 and is estimated by the inverse of the gap between the first 
two eigenvalues of H[y], Therefore, thanks to (|II40J) . we get that ||m||l 2 (r + ) < C\ T \ for 
some constant C > 0. Plugging this estimate together with (|II.40|) in (|II.38|) . we get (|II.39jl .D 

In the next proposition we have a regularity result with respect to the two variables (7, £). 
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Proposition II. 7 The functions (7, £) t— > and (7, £) i— > <p 7 £ are of class C°° in 

M. 2 . Moreover, we have : 

V 1} (7,O = l^(0)| 2 . (H.41) 



Using Proposition III. 51 and Remark 111.61 Proposition III. 71 is sufficient for achieving the 
proof of Proposition 111.41 

Proof of Proposition 111.71 In order to reduce the problem to a problem of an operator 
with a fixed domain, we define the bounded operator V[j] on L 2 (M + ) by : 

V[j]u = e 7xW *M, Vu G L 2 (M+). 

We define then the operator H[y,^] by : 

D(H[ 1 ,^]) = {ueB 2 (R + );u / (0)=0}, 
%,e] = V[-7]#[ 7 ,e]V[ 7 ]. 

Note that the domain of H [ 7 , £] is independent of 7 and £. Note also that H [7, £] is not 
self-adjoint but it has the same spectrum as H[y, £]. A fundamental state of H[y, £] is given 
by : 

We denote by if* ^ the orthogonal projector on Let us consider a point (70, £o)- We 

define the operator M : D{H[y,€\) xC — ► L 2 {R + ) x C by : 



M : 



#[7o,£o] - Ho V?o 



where /i = /i (1) (7o,£o) and v?o = <^ 7o ,£ Q . 

The operator M is invertible and its inverse Rq is given by 







£0 # + 
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where the coefficients of Rq are : 



E 



V[- l0 )R[ lo ,Z o )V[ lQ ) 



(11.42) 




V[-7o]<Po, 



(11.43) 



E, 



o 



(11.44) 




0. 



(11.45) 



The operator R[jo, £o] is the regularized resolvent which is equal to on M. ■ <p and to 
(#[70,6] -/^o) _1 on <pfr. 

Now we define, in a neighborhood of (7o,6o,A*q) the operator M(j, £, /i) by : 



The operator M(j, £, /i) is also invertible in a neighborhood of (70, £0, /io) and we denote its 
inverse by : 



• The coefficients of i?(7,£,/i) are C°° in a neighborhood of (7o,£o,A*o)- 

• A number fi is an eigenvalue of #[7, £] if and only if E~* (7, £, /i) = 0. 

Moreover, in a neighborhood of (70, £0, /-to), if A 4 i s an eigenvalue of H[j,^\, then 
V [y]E + ('y, £, fi) is a corresponding eigenfunction. 

Thus, in a neighborhood of (70,6), the eigenvalues of the operator H [7, £] are given by 
the solutions of the equation E + ~(j,£, n) = 0. By viewing the operator M(7,£,/i) as a 
perturbation of M , we can calculate the coefficients of _R(7,£,/i) and we obtain that : 



As the function E + (7, £, /1) is of class C°°, we can apply the implicit function theorem and 
get the existence of a number rj > and a function fi of class C°° such that : 





d^E + (70,6, /i ) = 1. 



V(7, e]7o - 77, 70 + ?/[x]£o - V, 6 + 77[, V/i g]/i ~V,Va + vl 



E + -( 7 ,e,/i) = 0^/i = M7,e)- 
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This proves that the functions (7,0 i— > /^(t, £) an d (7)0 l— > V?7,£ are °f class C 1 



□ 



C. Asymptotic behavior 

The asymptotic behavior at ±00 of the eigenvalue 0(7) with respect to the parameter 7 
is given in the following proposition. 

Proposition II. 8 There exist constants C ,7o > such that the eigenvalue 6(7) satisfies : 

1 - C o7 exp(-7 2 ) < 6(7) < 1, V 7 e [70, +oo[, (11.46) 



-7 2 < 0(7) < -7 2 + V 7 G]-oo,0[. (11.47) 



and 



Proof. We prove the estimate (jll.46|) . Note that by the min-max principle and Theorem lII.2l 

we get for any 7 > : 

^ (1) (O,O7))<0(7)<1- (H.48) 

The following estimate for the Neumann problem is obtained by Bolley-Helffepii (formula 
(A.18)) : 

|/x (1) ((U)-l| <C£exp-£ 2 , V£6[A+oo[, 
where C, A > are constants independent of £. Recalling (|II.15|) . the last estimate gives : 

k (1) (0,£(7)) - 1| < Co7exp- 7 2 , V 7 G [ 7 o,+oo[, 
where 70 = max(vC4, 1) and Co = 2C. Upon substitution in (jll.48)) . we arrive at the estimate 

(una. 

The relation (jII.15|) gives the lower bound 6(7) > — 7 2 . To get the upper bound in (|II.47|) . 
we use the function e 7 * (with 7 < 0) as a trial function for the quadratic form defining 
if [7,0], this which gives, 

H e Hl 2 (r + ) ^7 
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Therefore, we get by the min-max principle that /i^ 1 - ) (7,0) < — 7 2 H -. Recalling (111.61) . 

we get the upper bound in (jII47|) . □ 



D. Exponential decay of the ground state 



Using Agmon's technique (cf. Ref. Il5f ). we get the following decay result for the eigen- 
function ip T 

Proposition II. 9 For each e G]0, 1[ there is a positive constant C e such that, for all 7 G 1R, 
we have the following estimate for the eigenfunction tp 1 : 



exp ( e ) <y9 7 



<C £ (l + 7 - + 7- 



(11.49) 



ffi({*6K+;(*-5(7))>c e }) 
where we use the notation 7_ = max(— 7,0). 



Proof. Let us consider a function $ G if 1 (M + ). Given an integer iV G N, an integration by 
parts gives the following identity : 



N 



(eVr)' +|(*-e(7))eV 7 | 



alt + 7 |e $(0 V 7 (0) T - ^(iV)e 2 *^V 7 (^) 



2*(AT) 



@(7) ||e y^llz^o^]) + ||* e V ? 7lL2([ 0)A r ] ) 



(11.50) 



Let us recall that the eigenfunction </? 7 is strictly positive. It results then from the eigenvalue 
equation satisfied by y? 7 : 



V"(t) = {(t - £(7)) 2 - 6(7)) ^(t) > 0, V* g] v^t) + £(7), +oo[. 



Therefore, the function </?' is increasing on ]a/B(7) + £(7),+oo[. On the other hand, as 
V? 7 G H 2 {R + ), the Sobolev imbedding theore 
the monotonicity of ip' , we get finally that : 



<y9 7 G if 2 (lR + ), the Sobolev imbedding theorem gives lim <p'(t) = 0. Thus, combining with 



■;(*)< 0, VtG]v/e(7T + e(7),+oo[. 
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Taking iV > ^6(7) + £(7) and recalling that 6(7) < 1, the identity ()II.50[) yields the 
estimate : 

*n r 



(eVr + |(*-£(7))eV 2 



rft + 7 |e $(0 V 7 (0) 



- Il e ^t|Ii,2([ 0iJ v]) + || $ e V ? 7||l2([0,AT]) ' 



(11.51) 



To estimate the boundary term in (jll.51|) . we recall that ()11.38|) (with u — ip 1 and £ = £(7)) 



gives : 



K(0)| 2 <2 + ( 7 -) 2 . 
Therefore, the estimate (jll.51j) becomes : 



rN r 2 

/ (eV,)' +((^-e(7)) 2 -|$'r-l)|eV 2 
^0 



< 7- ^2 + ( 7 _) 2 e 2 * (0) . (11.52) 



Now we take $ as : 



We can then rewrite (|II.52|) as : 



$(t) = e 



(i-^(7)) S 



'*e[0,AT],(t-f(T))>ae 

where a f > satisfies : 



( e< V 7 )' 



1 w 

+ | e ¥>7 



tft < 7-^2 + (7-) 2 e e€(7) + e eae 



(11.53) 



a 2 — e 2 a f — 1 > 1. 



Notice that the first term on the right hand side of pi. 53)1 is effective only if 7 < 0. Coming 
back to the regularity of the function 6(7), the decay of 6(7) in (|II.47J) and the relation 
(jII.15|) . we get that the function £(7) is bounded for 7 < 0. Let us now take : 

C = sup £(7), C e = max (a e , e eCo ,e eclt ) . 

7<0 



The estimate ()H.53j) reads now as : 

(eV 7 )' 



1 I * I 2 

+ e ^7 



dt<a(i + 7- + (7-) 2 )- 



't&[o,N],(t-ai))>Cc 

Noticing that the above estimate is uniform with respect to N, we get (|II.49J) upon passing 
to the limit N -> +00. □ 
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Let us now recall that the regularized resolvent R{j] is the bounded operator defined on 
L 2 (R + ) by : 

{0 ; <p II (p~, 

(11.54) 
(i7[7]-e( 7 ))- 1 ; <f>±<p 7 , 

and extended by linearity Again, using the Agmon's technique, we get that this regularized 
resolvent is uniformly continuous in suitable weighted spaces. 

Proposition 11.10 For each S e]0, 1[ and rjo > 0, there exist positive constants Co, to such 
that, 

V 7 e [-770, m], Vn G L 2 (R + ; e^"*™)^) , u± ^, 

we /iave, 

e *(M(7))fl[ T ] u < Co||e 5 ^W)u|L 2m (11.55) 

III. PROOF OF THEOREM O 

In this section we prove Theorem 11.11 by comparing with the basic model introduced in 
the preceding section. We introduce a coordinate system (s, t) near the boundary dfl where 
t measures the distance to dQ and s measures the distance in dQ (cf. Appendix IX)) . 

Proposition III.l (Upper bound) 

Under the hypothesis of Theorem \I.1\ there exist positive constants C and h such that, 
Wh e]0, ho], we have : 

/i (1) («, 7, h) < h& (/i a ~ 1/2 (7o + Ch 1 ' 2 )) + Ctf' 2 . (III.l) 

Proof. We start with the easy case when a < \ and 7 > 0. Notice that in this case, given 
a constant C > 0, formula 46)1 gives the existence of ho > such that : 

|6 (h a - 1/2 { lQ + Ch 1/2 )) - 1| < exp(-/i 2Q - 1 ), V/i G]0, /i Q ]. (III.2) 

By comparing with the Dirichlet realization, the min-max principle gives : 

^\a n ,h)<X^(h), 
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where X^(h) is the first eigenvalue of the Dirichlet realization (on fl) of — (/tV — iA) 2 . Using 
the following upper bound for A«(„) (cf.Ref.3): 

X w (h) < h + Ch 3/2 , V/l G]0,1], 

together with (pTT2|) . we get fpTTT^ . 

We suppose now that 70 < if a < ^. Consider a point x G dfl such that 7(20) = 7o- We 
suppose that x = in the coordinate system (s,t) near the boundary (cf. Appendix lA"j) . 
Using this coordinate system we construct a trial function Uh, a supported in the rectangle 
Kh =] — h 1 / 4 , /i 1 / 4 [x [0, to[ following the idea of Helffer-Morame^ and Bernoff-Sternbergi. 
Since Xq is a minimum of 7, Taylor's formula up to the first order gives the existence of 
positive constants Cx, h such that, 

Vhe)0,ho), It(s) -7o| < Cxh 1 ' 2 in]-^ 1 / 2 ,^ 1 /^. 
Thus, given a trial function u supported in Kh, we have the following estimate : 

<do(«) < 9£jn(«)> V MM°], (ni.3) 

where 70 = 7o + Cih 1 / 2 . So it is enough to work with q^AVi- 

We introduce rj = h " 1 ^ 2 ^ and we choose now the following trial function : 



u h ,a = a 1/2 exp -i-TTjT- v h, a , (HL4) 



. 60?) s 

where a(s,£) = 1 — tn T (s) and 

«M = ^ 3/ V, (fc" 1/2 *) X(t) x / (/T 1/4 s) . (HI.5) 

The function x is a cut-off equal to 1 in a compact interval [0,to/2] and the function / G 
G o°(] - i is chosen such that ll/IU 2 (R) = L 

Note that the decay of in Proposition III. 91 gives 2 ^ : For every 5 > and G N, there 
exist positive constants C^s and /io such that, 

t k \y n {t)\ 2 dt < C h>s h' Sk , Vhe]0,h o ]. (III.6) 
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We work with the choice of gauge given in Proposition IA.21 Using formula (jA.3|) , we can 
write : 



|Sn|/2 rto 
-|3fi|/2 JO 



\hd t v h , a \ 2 + a~ 2 (h x ^(rj)-t(l 
\v h , a (s,0)\ 2 ds 



t 



-KAS) v h 



dsdt 



/|9fi|/2 
\v h>a (s,0)\ 2 c 
-\dfl\/2 

/\on\/2 rt 
/ [\(d t a-V 2 )v h J 2 + 2a- l ' 2 (d t a- 1 l 2 )(d t v Ka )v Ka + a~ 2 \d s v Ka \ 2 ] adsdt. 
■\dil\/2 JO 

(III.7) 

Recalling the expression of Vh, a (cf. (fllL5p ). we can replace the function x by 1 getting 
an exponentially small error on the right hand side of (|IH.7|) . thanks to the decay of tp v in 
Proposition III. 91 After a change of variables and using the decay of <p v in (IIII.6|) . the leading 
order term on the right hand side of (jIII.7|) is equal to : 

h [|^(*)| 2 + \(t - i{r)))^\ 2 dt] dt + V \ip v (0f 

and the error is of order 0(h 3 ^ 2 ). Therefore, we get constants C, ho > such that : 

QkXnM - h&( v )\ < C0 2 , Vh e]0, ho]. 

Using formula ()A.4|) and the decay of Lp v (Proposition III. 9)) . we obtain that the L 2 norm of 
Uh, a is exponentially close to 1 as h — > 0. The application of the min-max principle permits 
now to prove □ 



Remark III. 2 In the regime a g]|, 1[, we have, thanks to Proposition MI. 51 : 

9 (h a - 1/2 l0 ) = e + 6M 3 h a ~ 1/2 + OQi 2 *- 1 ). 

Substituting the above expansion in the upper bound hlll.l}) . we get the following upper bound 
for the eigenvalue //W(a,7,/i) ; 

/i (1) (a,7, h) < hG + 6M 37o //* +1/2 + 0(/i inf(3/2 ' 2a) ). 

We shall prove that this upper bound is actually an asymptotic expansion of /r 1 ) (a, 7, h) as 
h tends (see Remark W . 1 1\) . 
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Proposition III. 3 (Lower bound) 

Under the hypothesis of Theorem U.ll there exist positive constants C,C f and h such that, 
V/i e]0, ho], we have : 



{1 \a,-f,h) > hQ (h a - 1/2 l0 (l + C'h 1/4 )) -Ch 5/4 . (111.8) 



Proof. We follow the technique of Ref. |10| and we localize by means of a partition of unity 
to compare with the model operators in M? and K x R + . Let us explain the heuristic idea. 
A partition of unity permits to estimate the quadratic form n locally in small subsets of 
Q. Near the boundary, we obtain after a transformation of coordinates that the expression 
of q is to leading order asymptotics as that of the half-plane model. In the interior of 
Q, the expression of the quadratic form is actually like that of the entire-plane model. 
Let us introduce a partition of unity (xj) °f ^ 2 that satisfies 

l^'l 2 = 1; S I V ^'! 2 < +°°' su PPXj C D(zj, 1), 

3 3 

where for z6R 2 and r > 0, we denote by D(z, r) the disk of center z and radius r. 
We introduce now the scaled partition of unity : 

where eo and p are two positive numbers to be chosen suitably. Note that (x^) now satisfies : 

EWi a = 1 « ( IIL9 ) 

3 

J2\VXj\ 2 <Ce 2 h- 2 ?, (111.10) 

3 

supp xj C Q) := D(z$,e h<>), (III.ll) 
where C is a positive constant. We can also suppose that : 

either supp Xj H dtt = or z) G dQ. (III. 12) 
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Note that the alternative in (|III.12|) permits us to write the sum in (jIH.9|) under the form : 

£ = £+£• 



int bnd 



where the summation over "mt" means that the support of Xj do not meet the boundary 

while that over u bnd" means the converse. 

We have now the following decomposition formula : 

C>) = £CI(4")- ft2 £lll v 4MI 2 . (ni.i3) 

rV 

usually called the IMS formula (cf. Ref. 124( 1. We have now to bound from below each of 
the terms on the right hand side of fjlll. 13(1 . Note that (jHI.lOj) permits to estimate the 
contribution of the last term in (jlll.l3|l : 

h * Yl II I V *i HI' - Ce ° 2h2 ~ 2p \\ u \\ 2 , VueH^Sl). (111.14) 

3 

If Xj is supported in fl, then we have : 

CJ(x-u)= f \{hV-iA) X h j u\ 2 dx. 
Jr. 2 

Since the lowest eigenvalue of the Schrodinger operator with constant magnetic field in M? 
is equal to h, we get : 

> h [ \x-u\ 2 dx, V« G H\Q). (111.15) 
Jn 

We have now to estimate ^X^u) when Xj meets the boundary. It is in this case that we 
see the effect of the boundary condition. Note that, by writing Qh'lniXj 11 ) m t ne boundary 
coordinates, thanks to Proposition IA.1( there exists a positive constant C\ independent of 
h and j such that : 

/ |(W - iA) X jU\ 2 dx > (1 - de^hP) [ \(hV - iA) X jU\ 2 dsdt, \/u e H\tt), (111.16) 
Jn " Jrxr+ 



where A is the vector field associated to A by (|A.2[) . 

By a gauge transformation, we get a new magnetic potential A new j satisfying : 

A n ew,j = A — V0j , 
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A newJ (w) - Af in (w) <C\w\ 2 , w = (s,t), (111.17) 

where A\ in := \{—t, s) is the linear magnetic potential and C > is a constant independent 
of h and j. 

Given 9 > and any function u of support in K x R + , we get by the Cauchy-Schwarz 
inequality, 



/ (W - iA 

new. j) v ' {A 



,i ui,jf Aj in )v dsdt 



< h 20 / |(W - zA neu;i >| 2 + / - 



Writing A new j = A\ in + (A new j — Aj in ) and using (jIII.17|) . we get a positive constant C 
independent of h and j such that : 



|(/iV — iA new j)v \ 2 dsdt > (1 — h 



20\ 



- iAjjv 



dsdt-Ch- 2e II M 2 xN| 



(111.18) 



Let us recall that y^u is supported in the disk D^z 1 }, eoh p ). Upon noticing that 



|(W - iA)x h jU\ 2 dsdt 



(hV - iA newJ ) exp I -i^- x)^ 



h 



dsdt, 



we get by combining (|III.18|) (with v = exp ( — i-fi- ) together with (jIII.16)) . a constant 

C 2 > such that : 



|(W-zA)x-«| 2 c£r > (1 - C 2 e h p - C 2 h 2 



(hV-iAiJexp[-i^-} x)u 



-C 2 e h^- 29 



\XjU\ dsdt. 

Notice also, (possibly changing C 2 ), we have in D(zj,eoh p ), 



dsdt 
(111.19) 



TO*) > 7(4) " C26o^- 



Then, by putting, 



7j 



7 (4) - C 2 e h p 
I - C 2 h™ - C 2 t hP' 
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the estimate (jIII19|) reads finally : 



Lh 



-C 2 4h"- 2 »\\x h 1 uf. (111.20) 

Note that this permits to compare with the half-plane model operator and to get finally the 
energy estimate (cf. ([l.lOp ) : 



Cl,n(x» > {(1 - C 2 eoh> - C 2 h 2d )hQ(h^%) - Q&ti^"} \\xHl\ny (m.21) 
We substitute now the estimates (UlTSl) . (1111.14(1 . (1111.151) in 1)111. 13|) and get finally : 
CXnH>hJ2 [ \x*u\ 2 dx + hJ2®{h a - 1/2 l 3 ) [ \x]u\ 2 dx 

int ^ n bnd ^ U 

_C (h^- 2e + efh 2 - 2 " + h 1+ " + /i 1+2e ) ||u|| 2 , Vu G tf 1 ^). (111.22) 
As 7o is the minimum of 7, we can replace (|III.20|) by the estimate : 



-C 2e 2ft 4 "- 2e ||x>ll 2 , (HI.23) 

where 70 is defined by : 



7o 

7o 



' 1 - C 2 h 2e - C 2 e hP' 
We then get instead of (1111.221) : 



Cv»>^E / \x-u\ 2 dx + he (h a - 1/2 %)J2 j \Xju\ 2 dx 

int ^ n bnd ^ n 

_C {h Ap - 2e + tfh 2 - 2p + h 1+p + h 1+2e ) || M || 2 , Vm G F 1 ^). (111.24) 



The advantage of (jlll.22|) is that it gives a lower bound of the quadratic form n in 
terms of a potential, see however Section IIV1 

We choose now eo = 1, and we optimize by taking 2 — 2p = 1 + p = 4p — 29 (i.e. p = 3/8 
and 9 = 1/8) in (|IH.24|) . We obtain then ()HI.8|) by applying the min-max principle. □ 
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Proof of Theorem II. 11 The proof follows in principle from Propositions IIII.ll and IIII.3I 

Actually, in the regime a < |, we use further Proposition 111.81 while in the regime a > |, 
we use the continuity of the function 6(7) (Proposition III.4J) . □ 



IV. LOCALIZATION OF THE GROUND STATE 



We work in this section under the hypotheses of Theorem 11.31 Due to Theorem 11.11 we 
have in this case that 

lk&lfi<l. (IV.l) 

h->o h 

Then this gives, by following the same lines of the proof of Theorem 6.3 in Ref. ICt, the 
following proposition. 



Theorem IV.l Under the hypotheses of Theorem U . 31 there exist positive constants 5, C, h 
such that, for all h g]0, ho], a ground state u aa ^ of the operator P^Jq satisfies : 



Sd(x, dn) 

<'Xp ( -r-p j U a>1>h 



< C\\u 



L 2 (n) 



a,7,/i||L 2 (f2) ? 



and 



exp 



sd{x, an) 



hp 



<r- r~'/ 1 -min(l/2,/3) ||„, II 



where (3=1 — a if j < and a < ~, and [3=1/2 otherwise. 



(IV.2) 



(IV.3) 



Proof. Integrating by parts, we get for any Lipschitz function $ : 



QhXa ( ex P ( p ) u »n,h ) = /^ (1) (a ; 7, h) 



-h 



2-2/3 



|V$| exp ( —s ) u a ^ h 



2 

L 2 (n) 



(IV.4) 
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Let u = exp (A-) u a ,-y,h- Using the lower bound for q^'Aniu) in (|HL24|) together with the 
upper bound for fJ,^(a, 7, h) in we get from (jlV.4j) : 

V / (1 - 0(/i a - 1/2 ( 7o + C7i 1/2 )) - Ctfh 1 -*" - CJ&- 29 - 1 

int J V 

_ Ch min(p,20) _ /l l-2/9| V$ |2j x 

fend 7 ^ 



+ /l min(4p- 2e -l,l-2p) + fe l-2^| V$ |2j x [^f^ 



We choose p = j3 so that each is supported in a disk of radius t^hP . We choose also 6* > 
such that 4p — 26 1 — 1 > and we define the function $ by : 

$(x) = <5 max 

where 5 is a positive constant to be chosen appropriately. Note that 1 — Q(h a ~ 1 ^ 2; y ) decays 
in the following way : 

3 C , h > s. t., V7i e]0,ho], 

1 - 0(/i Q - 1/2 7o ) > Coh 2 "- 1 if 70 < and a < -, 

1 - 0(/i a - 1/2 7o ) > C otherwise. 

Thus we can choose eo and 5 small enough, so that we get finally the following decay : 

2 



n 



dx < C I \u ai y ; h\ 2 dx. 



This actually permits to conclude (|IV.2|) and, thanks to (|IV.4j) . 



€ln ( ^P s -^^u a>7>h ) < chr^w 



exp 



5d(x, dQ) 



(IV.5) 



L 2 {n) 



For a function u G let t) be defined by means of boundary coordinates (s, £) and 

equal to the restriction of u in fi to (cf. Appendix |XJ). Notice that : 

POO 

\u(s,0)\ 2 = -2 / {^(x(*)u(s,0)}xW"(s>*)^, 
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where x is the same cut-off introduced in (jIH.5|) . Integrating the above identity with re- 
spect to the variable s then applying a Cauchy-Schwarz inequality, we get after a change of 
variables the following interpolation inequality : 



u 



L 2 (an) 



where C is a positive constant depending only on Q. 

Applying again a Cauchy-Schwarz inequality, the preceding estimate gives : 
\\(hV - iA)u\\ 2 L2in) < 2qZl Q (u) + Ch\\uf L2(n) , \/ueH\n). 



In particular, for u = exp ^ 



Ua i h) we g e t (jlQ), thanks to (HE! and flEU). □ 



We study now the decay near the boundary. Let us consider a number (3 > and a 
Lipschitz function $o defined in Q. The function $ an d the number /3 will be chosen later 
in an appropriate manner. Choosing p — |, 9 — g and eo large enough, the energy estimate 
(1111. 22|) together with the upper bound ()111.1|) give the existence of a positive constant C 
such that : 

2 



M x h u h 



> h J2 I i 1 ~ (^^o) " Ch 1 ^ - h'-^ol 2 ) exp ^ 

int ^ Q ^ ' 

+/i V / [(e (^-^(x)) - 6 (/i Q - 1/2 7o)) - ^ 1/4 - ^- 2/3 |v<i»o| 2 ] 



dx 



exp 



o \ .A. 



Xj "a,7,)i 



where 7o = 7o + Ch 1 ^ 2 . The function 7 is extended to a small boundary sheath by means of 
boundary coordinates in the following way : 



7(x) = 7(5(2;)), Vx G Q 



to- 



In the case a < ~ and 70 = 0, thanks to Proposition 111.81 the difference between 
9 (/i a_1/,2 7(a;)) and 9 (/i a_1/,2 7o) decays in the following way : 

\/e > 0, 3C e > 0, Vx e (7 - 7or 1 ([e, +oo[), 9 (/i - 1 / 2 ^)) - 9 {h^ 1 ' 2 ^) > C £ . 
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In the case a < 5 and 70 < 0, we have a stronger decay : 

Me > 0, 3C £ > 0, Wx G (7 - 7o) _1 ([e, +oo[), 

e (/i a ~ 1/2 7(x)) - e (/i a " 1/2 7o) > c e h}- 2a . 

So by taking $0 i n the form : 

$o(x) = 5x(dist(x, <9f2))dist (a;, {x G <9f2; 7(2;) = 70}) , 

with 5 an appropriate positive constant and x is the same as in (|III.5J) . we get for each e > 
the following decay near the boundary : 

2 

dx <C e exp-^\\u a ^ h \\ 2 , V/iG]0A], (IV.6) 

with [3 = 1 — a if 70 < and a < |, and (3 = \ otherwise. This gives finally the decay in 
Theorem 11.31 

For the critical case a = | and 70 arbitrary, we define the function <3> by : 

$ (x) = 5 x(dist(x, dtt)) dist agm (x, {x G <9fi; 7(2:) = 70}) , 

where dist agm is the Agmon distance associated to the metric (6(7(0;)) — 6(70))+. We ob- 
tain then a similar decay result to (jlV.fiJ) . 

In the case when a > |, we need a finer energy estimate than (jIII.22)) . see however Re- 
mark |v^2 

V. TWO-TERM ASYMPTOTICS 

In this section we suppose in addition to the hypotheses of Theorem II. II that a > |. We 
give two-term asymptotic expansions for the ground state energy showing the influence of 
the scalar curvature and we finish the proofs of the remaining theorems announced in the 
introduction. 



dist(x, dfl)<tn 



$0 

6Xp ^ Ua,^,h 
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A. Upper bound 



We construct a trial function defined by means of boundary coordinates (s,t) near a 
point z G dfl. We suppose that z = in the coordinate system (s, t) and we denote by 
K o = K r(0), a = 1 — t^o an d ^(^o) = h a ~ 1 ^ 2, y(zo). We then define the trial function : 

■ £0?Oo)K 



M ft = exp -i 



h 1 / 2 



Vh(s,t), 



with 



v h (s,t) = h^ ie a 1/2 (t)v v{zo) (h- l 'H) X (t) ■ f {h-^s) 
and where the functions \ an d / are as m H 



(V.l) 



(V.2) 



We continue now to work in the spirit of Ref. 



L5j). 



lCl We work with the gauge given in Proposi- 



(V.3) 



tion IA.21 An explicit calculation, thanks to the decay of ip v (z ) fProposition 111.9(1 . gives the 
following lemma. 

Lemma V.l With the above notations, for each a G [|, 1] and 7 G C°°(<9f2; R) ; i/iere exist 
positive constants C, h such that, V/i G]0, /i ]j have the following estimate : 

- [ H h (C/V,(, )) x PA>„(*)) dt < Ch 13 
where the operators H h and U h are defined respectively by : 

H h = af (t (l - t^) - hV 2 av(zo))) 2 - h 2 a^d t (a d t 
{U h g){t) = h-^g(h-^t), Wg G L 2 (R + ). 
Proof. Note that in the support of Uh we have^i 

a = a + 0(/i 5/8 ), A 1 = -t(l- ^0) + 0(h g / s ) 

Then, thanks to formula (|A.3(1 (also cf. (|III.7(1 ) and the decay of (p v (z ) (Proposition |TT 
we get modulo 0(h 13 ^ 8 ) : 



"■() { \hd t v h \ + a 
+hV 2 V (z ) [ \v h (s,0)\ 2 . 



/(l-^o) -/r' \(//Co))) C, 



dsdt 



(V.4) 
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Integrating with respect to s, the right hand side above is equal to : 
h- 1/2 I a L 2 \d t (a - 1/ V,(,o)(^ V2 *)x(*)) 



+a~ n 



t ( 1 - -K ) - h^Zivizo)) ) V vM (h- 1/2 t) X (t) 



dsdt + hr)(z )\if v{zo) 



We can replace the function x hi the above expression by 1 getting an exponentially small 
error, thanks to Proposition 111.91 Thus, modulo a small exponential error, we rewrite the 
above expression as : 

2 



h 2 a \d t (U h <p v(zo) )\ 2 + a 2 



+h^ 2 V (z )\(U h ^ (zo) ) (0)| 2 . 



' I 1 - 5*0 ) - h^av^o)) ) {U\ v{zo) ) 



dsdt 



Notice that we have the boundary condition (U h (p v ( zo ^)'(0) = h~ 1 ^ 2 r](z )(U h cp ri ^ ZQ ))(0) . There- 
fore, integrating by parts, the above expression is equal to J R+ H h (U h <p v ( Zo ' ) ) x \ U h ^p ri ^ Zo ^) dt. 
Upon substituting in (|V.3jl . this finishes the proof of the lemma. □ 

Similar computations give also the following lemma. 



Lemma V.2 Under the hypotheses of Lemma W. 1\ there exist positive constants C,ho such 
that, Wh e]0,/i o ], we have : 

\\(H h - H h Q - H?)U h ^ M \\ L2{R+) < Ch 2 , (V.5) 

where the operators Hq and iff are defined respectively by : 
HZ = -h 2 d? + (t-hV 2 t(r)(z ))) 2 , 

E\ = 2tn (t - h^faizo))) 2 - K t 2 (t - hWfaizo))) + h 2 K d t . 
Let us denote by (cf. $FE$) and (|TT25jl ) : 

M 3 (~,l(zo)\ = ^3(7(^0)), M 3 (a,j(z )) = M 3 for a > i. 
The next lemma permits to conclude an upper bound for the eigenvalue fj,^(a, 7, h). 
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Lemma V.3 Under the above notations, there exist positive constants C, ho such that, when 
h e]0, ho], we have the following estimate : 



fioW " {©fate))) " 2M 3 (a, 1 (z o ))K O 2 }\\u h \\l 2(n) 



where e a = inf (13/8, 2a + ~) /or a > | and ex/ 2 = 13/8. 



< C7i e <*, 



Proof. Notice that in the support of Uh we have 7(2) = 7(2:0) + (^(/i 1 ^ 8 )- Then this gives 



%?A^u h )-q^S\u h ) = 0{h^). 



In view of Lemmas IV . II and IV .21 we get the following estimate : 



ClnK) " / {Hi + #o) {U h V vM ) x U\ v{zo) dt 



< Ch 13 ' 8 . 



(V.6) 



We note also that we have the following relations : 

(U h )*H*U h = h {-d? + (t - ^ ))) 2 } 



(V.7) 
(V.8) 



where the operator Hi is defined by : 



H l = {t- £(r/U)))) 3 - £(^o)) 2 (t - av(zo))) + ft. 



By defining K 3 (a, h) := f R Hip^^ ■ (p v ( zo )dt, the estimate (|V.6J) reads as : 
KlnM - {hQ^zo)) + IU(a,h)KohV 2 }\ < Ch 13 / 8 . 



(V.9) 



Now, for a — |, we get by using (|II.3J) that K 3 {\, h) = — 2M 3 (|, 7(^0)). For a > |, thanks 
to Propositions 111.41 and 111.51 we get that 



2a-l\ 



K 3 (a,h) = -2M 3 + 0(h 

Finally, the decay of <p n ( Zo ) m Proposition 111.91 gives that H^Hl 2 ^) is exponentially close to 
1. This achieves the proof of the lemma. □ 
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The min-max principle gives now, thanks to Lemma IV. 31 an upper bound for //W(a, 7, h). 
Under the hypothesis of Theorem II. 4[ we take z such that 



(«r - 37) (z ) = (Kr ~ 37) max 

and we use the expansion (cf. (|II.27|) ) : 

Q(v{z )) = 6 + 6M 37 (^o)/i 1/2 + 0{h). 
Therefore, (|V.9|) gives the following upper bound : 

//«(1,7, /i) < /^e - 2M 3 (k t - 3 7 ) max /i 3/2 + 0(/i 13/8 ). (V.10) 
Under the hypothesis of Theorem II .5| we choose z such that n r (z ) = (« r ) ma3c . 

B. Lower bound 

As in the proof of Proposition IIII.3I we consider a standard scaled partition of unity2i 

(XjViVs) - eZ 2 of 1R 2 that satisfies : 

E wi 2 = 1, E i v w(*)i 2 ^ c ^ 1/3 ' ( v - n ) 

supp Xijfcl/8 C j/i 1/6 + [-/i 1/6 , /i 1/6 ] 2 - (V.12) 
We define the following set of indices : 

J l{h) '■= {3 e Z 2 ; BuppXj, h i/6 n ft ^ 0, dist(supp x ijfc i/a, 5ft) < r(/i)}, 

where the number r(h) is defined by : 

T {h) = h 5 , with - < 5 < -, (V.13) 
6 2 

and the number 5 will be chosen in a suitable manner. 
We consider also another scaled partition of unity in R : 

<r(A)(*) + (*) = X > IV4, T (A)(*)I <^y J = 0, 1, (V.14) 
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SUpp^ , r (h) C [-^-, + Oo[, SUpp^i iT ( fc ) Cj - 00, 



(V.15) 



Note that, for each j 6 ^r(h)> ^he function V , i, r (h)(^)Xj,/ii/6('S, £) could be interpreted, by 
means of boundary coordinates, as a function in f2. Moreover, each V , i,r(fi)(£)Xj,/i 1 / 6 ( s J * s 
supported in a rectangle 

near dQ. The role of 5 is then to control the size of the width of each rectangle K(j,h). 
Due to the exponential decay of a ground state away from the boundary ( Theorem IIV.1|) . 
we get the following lemma. 

Lemma V.4 Suppose that a > \. With the above notations, a L 2 -normalized ground state 
u a,-y,h of the operator Pff'J n satisfies : 



< Ch 5/3 . 



(V.16) 



The proof of ( V.16) follows the same lines of that in Ref . \\\ (Formulas (10.4), (10.5) and 
(10.6)). 

For each j £ J\h)i we define a unique point Zj G dfl by the relation s(zj) = Sj. We denote 
then by Kj = K r (zj), cij(t) = 1 — Kjt, A^(t) — — t (l — fftj), and = j(zj). 
We consider now the fc-family of one dimensional differential operators : 



H hJtk = -tfa-'dtiajdt) + (1 + 2 Kj t){hk - A j ) 



(V.17) 



where A; is a real parameter. We denote by H^'J'® the self-adjoint realization on 
L 2 (]0, h s [] aj(t)dt) of Hhj,k whose domain is given by : 



D{Hlik) = {« e H 2 (]0,h d [);v'(0) = h a j jV (0), v{h 6 ) = 0}. 



(V.18) 



The parameter jj is defined by : 
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where e(h) = if the function 7 is constant; if 7 is not constant, then there are constants 
C, ho > such that : 

\e(h)\ < Ch l '\ V/te]0,/i o ], 



We now introduce : 



/4(«,7,/0 := M infSp(^f fe ). (V.19) 



fcGM 

We have now the following lemma. 

Lemma V.5 For each a G [h, +00 [, we /iawe under the above notations : 



7, *») > ( inf Ai}(a,7,/») ) + 0{h^). (V.20) 

J &J T(h) 



Again the proof follows the same lines of Ref. |lOj (Section 11), but let us explain briefly 
the main steps. We express each term o[i})i,T(h)Xj,h' l / <iU a,'Y,h) in boundary coordinates. 
We work with the local choice of gauge given in Proposition IA.21 We expand now all 
terms by Taylor's formula near (sj,0). After controlling the remainder terms, thanks to 
the exponential decay of the ground states away from the boundary, we apply a partial 
Fourier transformation in the tangential variable s and we get finally the result of the lemma. 

We have now to find, uniformly over k G R, a lower bound for the first eigenvalue 



Hi(k; a, 7, h) of the operator H^'J' k . Putting (3 = k,j, £ = —h l l 2 k and 77 = 7,-, we get by a 
scaling argument : 

Hx(k; a, 7, h) = hfii(H%jjg), 
where fj,i(H^g) is the first eigenvalue of the one dimensional operator : 

Kjf = ~9t +(t- 2 + Ph 1/2 (1 - ph x IH)- x dt 

+2l3h 1 'h( K t-i-l3h 1 l 2t ^j -ph^H 2 (t-0+P 2 hj, (V.21) 

whose domain is defined by : 

D(Kpf) = {« e ff 3 ao ) /i'- 1 / 2 D;u , (o) = fc^VCo), u(h 5 ^ 2 ) = 0}. 
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We have then to find (when tj,/3 g] — M,M[ and M a given positive constant), uniformly 
with respect to £ G M, a lower bound for the eigenvalue /^(if^'^P). The min-max principle 
gives the following preliminary localization of the spectrum of the operator H^'^ : 

Lemma V.6 For each M > and a G [|,+oo[, i/iere exist positive constants C,ho such 
that, 

Vr],{3 e}- M,M[, V£GM, \/he]0,h ], 

we have, 

HiKS°) - NiKt) I < (l + N {Kl D )) , (V.22) 

where, for an operator T having a compact resolvent, Pj{T) denotes the increasing sequence 
of eigenvalues ofT. 

Remark V.7 Note that the min-max principle gives now that 

where, for fj G R, fi^'(fj,^) is the increasing sequence of eigenvalues of the operator H[r),£] 
introduced in Ul.fy) . 

The following lemma deals with the case when £ is not localized very close to i(h a - l / 2 r]). 



Lemma V.8 Suppose that 5 G]l/4, l/2[. For each a > \, there exists p G]0,<5 — |] ; and for 
each M > 0, there exist positive constants (, h > such that, 

Vt7,/3 g]-M,M[, such that \i-i(h a ~ l ' 2 r])\>Ch p , Vhe}0,h o ], 

we have, 

Vi(Kjf) > e(^- 1/2 r/) + (V.23) 

Proof. It is sufficient to obtain (jV.23)l for p,^{h a ~ l / 2 r], £), thanks to Lemma IV.6I and 
Remark IV.71 We start with the case when a = \ and r\ g] — M, M[. Writing Taylors 
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formula up to the second order for the function £ i— > ^'(rj, £), we get positive constants 
6>, Ci such that when |£ — £(r})\ < 9, we have : 

Then by taking ( such that Ci£ > Co, where Co > 1 is a constant to be chosen appropriately, 
we get when (h p < |£ — £(77) | < 0, 

^ (1) (^0>©(^) + Co^, 

where p is also a positive constant to be chosen later. When |£ — £(77) | > we get a positive 
constant eg such that : 

Then by choosing h such that (0^0 < e f' we S e ^ f° r l£ ~~ ^l 7 ?)! — (h p an d ^ ^0] : 

» {1) (v,0>®(v)+(oh 2p . (V.24) 

We treat now the case when a > 1/2. Note that the min-max principle gives uniformly for 
all £ e K and 7) e] - M,M[, 

^(1)^-1/2^ > (1 _ ^^-1/2)^(1) ( Q) £). 

Then using ()V.24|) for 77 = and p = inf (5 — 7, a — |), we can choose Co large enough so 
that we have for |£ — £ | > C^, p : 



<2 
2 

To finish the proof, we replace £0 by £(/i a-1 / 2 7/) getting an error of order 0{h a ~ l l 2 ). □ 



Now we deal with the case when |£ — £(/i Q 1 ^ 2 7/)| < (h p . Let i] = h a 1 ' 2 r). We look for a 
formal solution (//, f^g^) of the spectral problem 

tf^VX = /4Sk (o) = h a - 1/2 fZ;ls(o), (V.25) 



43 



in the form : 



n = d + dt (e - m) + d 2 (e - ^)) 2 + 4^ 1/2 , 
/£L = u ° + (t- ui + - am 2 u 2 + /i i/2 n 3 , 



(V.26) 
(V.27) 



where the coefficients do,di,d 2 , d 3 and the functions uo, ui, u 2 , w 3 are to be determined. We 
expand the operator if^'^P in powers of (£ — £(f})) and then we identify the coefficients of 
the terms of orders (£ — (j = 0, 1, 2) and h 1 ^ 2 . We then obtain for the coefficients : 

d = Q(fj), u = tp.^ 
d 1= 0, Ul = 2R{r}}{(t-av))Vv} 
d 2 =: d 2 (a,r]) = 1 - 2 J u+ (t - ^(fj))(pfjUidt, 

(V.28) 

u 2 = R[fj] {4 (t - £(fj)) R[fj] [(t - <p n ] - d 2 ] 

d 3 =: d 3 (a, 77) = (3 J R+ iff, {d t + (t - £(r?)) 3 } ffjdt, 

k % = [/? (3t + (t - £(??)) 3 - £(r/) 2 (f - m)) ~ M uo- 

Using the function x( ^-1/2 ) fh'p g (where x is the same as in ()IIL5p ) as a quasi-mode, we get 
by the spectral theorem, thanks to the decay results in Propositions 111.91 and III. 101 and to 
the localization of the spectrum in Lemma IV. 61 the following lemma. 

Lemma V.9 Suppose that 5 ~[. For each M > and a G [|, 1], there exist positive 
constants C > 0, h such that, 



we have, 



V77, /3e]-M,M[, V£ such that |f - £(77) | < (h p , V7i e]0, fr ] 



(V.29) 



< C [^1^-^)1 + ^/2] _ 
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where d2{ot,rj) and d^{a,rj) are defined by \ V. 28^1 respectively. 

Hence we have obtained by this analysis a lower bound for the first eigenvalue //-^(a, 7, h). 
We complete the picture by showing that the term d 2 (a,rj) is positive. 

Lemma V.10 For each a G [~, +oo[ and M > 0, there exists a positive constant ho such 
that : 

d 2 (a,r])>0, Vhe]0,h }, Vr/G]-M,M[. 

Proof. It is actually sufficient to prove the conclusion of the lemma when a = \. If a > |, 
we replace d 2 (a,rj) by its approximation up to the first order, thanks to Proposition 111.41 
and we obtain that : 

d 2 (a, V ) = d 2 (± 0) + O(h a - 1 / 2 ), 

which gives the lemma. 

For the particular case a = ~, we show that : 

which is strictly positive. □ 

We are now able to conclude the asymptotics given in Theorems 11.41 and 11.51 First we 
choose S — JL When a > | we replace B(^) and d^ by their approximations up to the 
second and first orders respectively, thanks to Propositions III. 51 and 111.41 For a = | we get 
by (jll.21j) that d 3 is indeed equal to — 2M 3 (|,7). 

Remark V.ll When a g]|, 1[ and when the function 7 is not constant, we get from the 
above analysis that the upper bound in Remark \lll.2\ is indeed an asymptotic expansion, and 
we achieve therefore the proof of Theorem \I.£\ 

We get also that the quadratic form q°^\ n can be bounded from below by means of a potential 
W : 

qt'XM > I W{x)\u{x)\ 2 dx, Wu e H\Q), 
Jn 
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where W is defined for some positive constant Co by : 

{ h -if dist(x,9fi) > /i 1/6 

W(x) = { 

{ h& + QM 3l (x)h a+1 / 2 - C /i inf(3/2 - 2a) ; if dist(z, dQ) < h 1 ^. 
Then, as in Section \TV[ we get by Agmon's technique that a ground state decays exponentially 
away from the boundary points where 7 is minimum and hence we have completed the proof 
of Theorem \I.3[ 

Remark V.12 Note also that the above analysis permits, under the hypotheses of Theo- 
rems \1.4\ and \l.d\ to bound the quadratic form q°^\ n from below using a potential W defined 
either by means of the function k t — 37 ( when a = 1) or by the scalar curvature k t ( when 
7 is constant). Then, by using Agmon's technique, we finish the proofs of Theorems \I.J\ and 




VI. CONCLUSION 



1G has allowed us to understand the role of 



The systematic analysis in the spirit of Ref. 
the boundary condition imposed by De Gennes. We have extended in Theorems 11.41 and 11.51 
the expansion announced by Pan2£ in the particular case when a = 1 and 7 is a positive 
constant. However, there is a specific difficulty when 7 is negative. We have not been able 
to obtain the localization of the ground state when a < 1/2 and 70 > 0. This is strongly 
related to the question of the localization of the ground state of the Dirichlet realization of 
the Schrodinger operator with constant magnetic field which is open. Finally, in the spirit 
of Refs. Il2ll26lj27l we hope to apply this analysis to the onset of superconductivity and to 



complete the analysis of Ref. 



23 (cf. Ref.llSD. 
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APPENDIX A: COORDINATES NEAR THE BOUNDARY 

We recall in this appendix well-known coordinates that straightens a portion of the bound- 
ary <9f2. Let sg]- ^p^, 4^-] I— > M(s) G <9f2 be a regular parametrization of cftl For each 
x G VL and e > we denote by : 

t(x) = dist(x, dQ) and Q e = {x G Q; dist(x, dfl) < e}. 

Then there exist a positive constant t > depending on Q such that, for each x G f2 to , we 
can define the coordinates (s(x),t(x)) by : 

t( x ) = \x-M(s(x))\, 

and such that the transformation : 

ij :Q tQ 3x^ (s(x),t(x)) G S| aQ | /27r x [0,t [ 

is a diffeomorphisim. The Jacobian of this coordinate transformation is given by : 

a(s,t) = det(Dtp) = 1 -tK T (s). (A.l) 

To a vector field A = (A 1 ,A 2 ) G C°°(n;lR 2 ), we associate the vector field A = (A 1 ,A 2 ) G 
C°°(§| (9n | /27r x [0,i [) by the following relation : 

Aids + A 2 dt = A-idxi + A 2 dx 2 . (A. 2) 

We get then the following change of variable formulas. 
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Proposition A.l Let u G H 1 ^) be supported in Q to - Then we have : 

I |(W - iA)u\ 2 dx = [ \\(hd t -iA 2 )v\ 2 + a- 2 \(hd s -iA 1 )v\ 2 \ adsdt. (A.3) 



and 



u{x)\ 2 dx = / \v{s,t)\ 2 adsdt, (A.4) 

where v(s, t) = t)). 
We have also the relation : 

(d Xl A 2 — d X2 Ai) dx\ A dx 2 = [d s A 2 — d t A^ a~ x ds A dt, 

which gives, 

curl A = (1 — tK r (s)) curl A 
We give in the next proposition a standard choice of gauge. 

Proposition A. 2 Consider a vector field A = (Ai,A 2 ) G C 00 ^;!! 2 ) such that curL4 = 1. 
For each point x$ G dQ, there exist a neighborhood V xo C £lt °f x o an d a smooth real-valued 
function <fr Xo such the vector field A new := A — V0 Xo satisfies : 

A l new = -t(l- ^k t (s)J and A 2 new = in V XQ , (A.5) 

with A new = (A\ ew , A 2 



\ 2 ) 

new J 
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Actually we shall need this decay only when a < 1/2 and 70 < 0. 

Actually, if zq is a point of maximum of k t , the remainder is better and of order 0(/i 3 / 4 ) for 
the first term. 

We take a partition of unity associated to squares instead of discs. 
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